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On  the  basis  of  general  assumptions  that  are  pertinent  to  penetration  modeling,  we  inves¬ 
tigated  the  effect  of  layering  of  metallic  shields  on  their  protective  characteristics  against 
sharp-nosed  rigid  impactors  and  established/explained  a  number  of  general  features  of 
the  penetration  phenomenon.  The  main  results  can  be  formulated  as  follows:  (i)  any  lay¬ 
ering  has  an  adverse  effect  on  protective  properties  of  the  shield  against  perforation,  (ii) 
increasing  the  number  of  layers  with  the  same  thicknesses  impairs  protective  properties 
of  the  shield.  These  assertions  are  validated  by  analyzing  results  of  experiments  and 
numerical  calculations  available  in  the  literature  that  allow  comparing  protective  proper¬ 
ties  of  monolithic  and  layered  shields. 

©  2010  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Effect  of  layering,  i.e„  replacing  the  monolithic  plate  by  several  plates  in  contact  with  the  same  total  thickness,  on  the 
ballistic  properties  of  metallic  shields  against  rigid  projectiles  has  attracted  considerable  interest  of  the  researchers  in  the 
field  for  a  long  time.  More  detailed  information  can  be  found  in  the  surveys  [1,2],  monograph  [3]  and  recent  publications 
[4-12], 

Quite  a  few  publications  examine  sharp-nosed  impactors  [4-10,13-18]  that  are  considered  in  this  study.  Most  of  these 
studies  report  impairing  of  protective  properties  of  the  shield  as  a  result  of  layering  which  is  observed  in  experiments 
and  numerical  calculations.  Using  a  semi-empirical  model  Radin  and  Goldsmith  [14]  proved  that  the  ballistic  limit  velocity 
(BLV)  of  a  two-layer  shield  with  layers  in  contact  and  having  the  same  thicknesses  is  smaller  than  the  BLV  of  a  monolithic 
shield  manufactured  from  the  same  material  and  having  the  same  total  thickness.  To  the  best  of  our  knowledge  the  latter 
study  is  the  only  one  which  employed  analytical  approach. 

In  present  study  we  investigate  the  effect  of  layering,  namely,  the  thicknesses  and  the  number  of  the  layers,  on  the  pro¬ 
tective  property  of  a  shield.  This  study  is  based  on  the  approach  that  we  proposed  and  used  recently  to  concrete  shields  [19]. 
The  special  feature  of  this  approach  is  that  it  does  not  require  a  closed  mathematical  model  for  calculating  BLV  of  monolithic 
and  layered  shields.  We  showed  that  many  characteristics  of  layered  shields  can  be  explained  and  predicted  on  the  basis  of 
general  properties  of  the  functions  which  determine  penetration  models  without  invoking  particular  expressions  for  these 
functions.  This  is  very  important  since  qualitative  behavior  of  the  characteristics  in  penetration  mechanics  is  generally  deter¬ 
mined  using  simple  models  which  include  approximations  of  the  experimental  data. 
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Nomenclature 

D  diameter  of  the  cylindrical  part  of  the  impactor 

g  function  defined  in  Eq.  (13) 

G  function  in  Eq.  (2) 

H<v/'),  H(v)  elements  of  the  Hessian 
fci .  k2  parameters  of  models 
m  mass  of  the  impactor 

N  number  of  plates  in  the  layered  shield 

T  thickness  of  a  shield 

Tmin  lower  bound  of  the  domain  of  function  G(T) 
Tsum  total  thickness  of  the  layers  in  the  shield 
T®  thickness  of  the  jth  plate 
vbt  ballistic  limit  velocity  of  the  monolithic  shield 
Vres  residual  velocity  for  the  monolithic  shield 
Vw  ballistic  limit  velocity  of  the  layered  shield 

V®  BLV  of  the  jth  plate 

V imp  impact  velocity  of  the  layered  shield 

l/P'p  impact  velocity  of  the  jth  plate 
Vres  residual  velocity  of  the  layered  shield 
V®  residual  velocity  of  the  jth  plate 
a,  f!  parameters  of  the  model  in  Eq.  (16) 
r] res,  f/w  parameters  in  Eq.  (17) 

<5  parameter  defined  in  Eq.  (7) 

0  nose  angle  for  the  conical-nosed  projectile 

A  criterion  of  optimization  in  Eq.  (B.l) 

/(0,  /(j  parameter  of  the  model  in  Eq.  (10) 

BLV  ballistic  limit  velocity 


2.  Mathematical  models  of  the  layered  shield 


Consider  a  high  speed  normal  penetration  of  a  rigid  sharp-nosed  striker  into  a  shield  consisting  of  N  metal  plates  in  con¬ 
tact.  The  plates  are  numbered  in  the  direction  of  the  impactor’ s  motion  as  1.2 . N  and  the  thickness  of  the  ith  plate  is  T®. 

Based  on  the  frequently  used  assumption  that  the  layers  do  not  interact  and  are  perforated  independently,  the  residual 
velocity  at  the  jth  plate  V®.  is  equal  to  the  impact  velocity  at  the  (j+1  jth  plate  Vj"  for  j  =  1  4-  N  -  1 .  In  addition,  the  impact 
velocity  of  the  shield  Vimp  is  equal  to  Vj®p,  and  the  residual  velocity  of  the  shield  Vres  is  equal  to  V®1  =  v!^p1)  (the  notation 
's  introduced  here  for  convenience).  Hereafter,  the  superscript  in  the  brackets  denotes  the  ordinal  number  of  the  layer. 

We  use  the  following  relationship  between  the  impact  velocity  V®p,  the  residual  velocity  V®  and  the  BLV  of  the  plates 

[20]: 


1/0)2  _  V(j)2  v(j)2 

imp  res  v  bl  ’ 


(i) 


where  the  BLV  is  defined  as  the  minimum  impact  velocity  of  the  impactor  required  for  perforating  a  shield  whereby  impactor 
emerges  from  the  shield  with  a  zero  residual  velocity. 

Since  the  plates  are  manufactured  from  the  same  material,  the  BLV  of  a  plate  is  a  function  of  its  thickness: 

V®2  =  G(T®),  G(0)  =  0,  j  =  1,2, ...  ,N,  (2) 

where  G  is  an  increasing  function. 

Summing  Eq.  (1)  from  j  =  1  to  j  =  N  yields: 


V2 

imp 


-v2 


]Tc(r®), 


j=i 


(3) 


and,  consequently,  the  expression  for  the  BLV  of  the  layered  shield  Vw  reads: 

V2i  =  £c(T®). 
j'=i 


(4) 


In  the  case  of  the  monolithic  shield  with  the  thickness 
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Tsum  =  Tm  +  T(2)  +  ■  •  ■  +  Tm,  (5) 

the  relationship  between  the  BLV,  vbh  and  the  thickness  of  the  shield  can  be  obtained  from  Eq.  (4)  substituting 
Vbi  =  Vu,  N  =  1  and  T(1)  =  Tsum: 

v2bl  =  G(Tsum)  =  c(j2T^  j.  (6) 


3.  Comparison  of  monolithic  and  layered  shields 

In  this  section,  we  compare  the  magnitudes  of  the  BLVs  of  the  layered  shield  and  of  the  monolithic  shield  having  the  same 
thickness.  Clearly,  one  can  consider  the  difference: 

<5  =  <  -  V2bl  =  G  ( jr  T®)  -  jr  G(T® ),  (7) 

\j'=i  /  j=i 

rather  than  the  difference  of  the  BLVs,  vM  -  Vb/.  The  main  goal  of  this  Section  is  to  prove  the  following  claim. 

Claim.  Let  G(T )  be  a  convex  (convex  downwards),  twice  differentiable  function  which  is  defined  in  the  interval 
T  3=  T min  3=  0  and  satisfies  the  inequality: 

TminG'(Tmin)  —  G(Tmin)  >  0,  (8) 

where  Tmin  is  given.  Then  5  >  0  for 

r0)>rmin,  t®  >  o,  j  1.2 . n.  (9) 

The  proof  of  above  claim  is  given  in  Appendix  A. 

Introducing  the  parameter  Tmin  allows  us  to  separate  the  range  of  the  considered  plate  thicknesses.  If  Tmin  =  0,  then 
G(0)  =  0  as  can  be  deduced  by  simple  physically  reasoning.  When  Tmin  >  0  function  G(T)  is  not  defined  for  T  <  Tmin  and  for¬ 
mal  extension  of  this  dependence  for  T  <  Tmin  does  not  have  physical  sense  and  often  implies  that  G(0)  ¥=  0.  Geometrically, 
Eq.  (8)  can  be  interpreted  as  follows.  The  tangent  at  the  point  T  =  Tmin  to  the  curve  determined  by  the  function  w  =  G(T )  in 
the  T,  w  plane  intersects  the  axis  T  =  0  at  the  point  w0  <  0.  In  particular,  this  condition  reduces  to  the  inequality  G(0)  ^  0  if 
Tmin  =  0.  Eq.  (8)  can  be  changed  by  continuing  the  function  w  =  G(T)  up  to  T  =  0  and  retaining  the  convexity  so  that  it  inter¬ 
sects  the  axis  T  =  0  at  w0  ^  0. 

If  G(T)  is  a  linear  function, 

G(T)  =  fitT  +  f.  l0,  (10) 

then 

5  =  ^  g r®  +  ^  -  £>,1®  +  fi0)  =  n0{\  -N).  (li) 

Thus  S  >  0  if  f.i0  <  0  and  5  =  0  if  fi0  =  0. 

Consequently,  the  monolithic  shield  is  superior  over  any  layered  shield  with  the  same  total  thickness  if  G(T)  is  convex 
function  and  conditions  given  by  Eqs.  (8)  and  (9)  are  satisfied  or  G(T)  is  a  linear  function  with  negative  constant  term,  /r0. 
If  G(T)  oc  T,  then  layering  does  not  effect  ballistic  properties  of  the  shield. 

4.  The  worst  layering  for  a  given  number  of  layers 

According  to  the  above  analysis,  ballistic  properties  of  the  shield  can  be  impaired  by  layering.  Let  us  estimate  the  max¬ 
imum  reduction  of  the  ballistic  performance  of  the  shield.  To  this  end,  we  will  show  that  this  occurs  when  the  layers  have  the 
same  thicknesses.  We  assume  here  that  the  number  of  the  plates  in  the  shield  N  and  the  total  thickness  of  the  shield  Tsum  are 
given.  The  problem  is  to  minimize  Vw  in  Eq.  (4)  under  constraints  given  by  Eqs.  (5)  and  (9).  We  also  assume  that 

Tsum/N  3*  T mm,  i.e.,  the  thicknesses  T®  =  Tsum/N,  j  =  1,2 . N  are  permissible. 

Investigation  of  this  optimization  problem  with  convex  function  G(T),  implies  (see  Appendix  B)  that  the  worst  layering 
occurs  when  the  monolithic  shield  is  replaced  by  the  identical  plates  having  the  same  total  thickness.  If  G(T )  is  a  linear  func¬ 
tion  given  by  Eq.  (10)  then 

K  =  X>,r®  +  Ho)  =  Hjsun,  +  H0N,  (12) 

J=1 

i.e.,  sizes  of  the  layers  do  not  affect  ballistic  properties  of  the  shield. 
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5.  Effect  of  the  number  of  layers 

Let  us  determine  the  dependence  of  the  BLV  on  the  number  of  the  plates  in  the  shield  assuming  that  all  plates  have  the 
same  thicknesses.  The  latter  assumption  implies  the  following  expression  for  the  BLV  of  the  layered  shield: 

i?bl  =  g(N)  =  NC(Tsum/N).  (13) 

Assuming  that  N  varies  continuously  we  can  calculate  the  derivative: 

g'(N)  =  C(T)  —  TC'(T),  T  =  Tsum/N.  (14) 

It  is  proved  in  Appendix  A  that  Eq.  (8)  implies  Eq.  (A.l)  when  G(T)  is  a  convex  function.  Equation  (A.l),  in  its  turn,  implies 
that  g'(N)  <  0  if  Tsum/N  >  Tmin.  Therefore  g(N)  decreases  when 

1  <  N  ^  Tsum/Tmin.  (15) 

The  conclusion  is  that  the  BLV  decreases  with  the  increase  of  the  number  of  layers  if  all  plates  have  the  same  thicknesses 
and  G(T)  is  a  convex  function. 

In  the  case  of  a  linear  function  G(T),  Eq.  (12)  implies  that  the  BLV  decreases  with  the  increase  of  the  number  of  layers 
(N  <  -iqTsum //(„),  independent  of  their  thicknesses  if  /.i0  <  0  and  does  not  depend  on  of  the  number  of  layers  if  fi0  =  0. 


6.  Discussion  of  assumptions  and  results 


There  are  only  a  few  publications  which  can  be  used  for  the  analysis  of  behavior  of  the  function  G(T). 
The  model  of  Ohte  et  al.  [21]  for  conical-nosed  impactors  implies  that 


v 


2 

bl 


k,T3 


a  =  1  +  2.9 


1 

a’ 


(16) 


where  m  and  D  are  the  mass  and  the  diameter  of  the  impactor,  correspondingly,  0  is  the  nose  angle  for  the  conical-nosed 
projectile,  T  is  the  thickness  of  the  shield,  and  the  coefficient  p  depends  on  the  system  of  units.  Therefore,  G(T)  =  /qT3  is  a 
convex  function. 

Gupta  et  al.  [22]  describes  the  dependence  between  the  BLV  of  the  ogive-nosed  impactor  and  the  thickness  of  the  alumi¬ 
num  plate  as  Vu  =  k2T089,  where  l<2  is  a  constant.  In  this  case,  G(T)  =  k2P  78  is  a  convex  function.  Functions  G(T)  in  Fig.  la, 
which  we  plotted  using  experimental  results  [6,14],  are  convex  and  satisfy  the  condition  given  by  Eq.  (8).  Processing  of  the 
experimental  results  [23]  (Fig.  lb)  yields  a  linear  function  G(T),  with  /,i0  <  0  in  three  out  of  four  cases  and  /t0  =  0  in  one  case. 

Results  of  the  experiment  and  numerical  calculations  [5-8,13-15]  used  in  the  analysis  of  our  theoretical  predictions  are 
summarized  in  Tables  1-8.  All  experiments/calculations  in  each  table  are  obtained  for  the  impactor  having  the  same  shape 
against  the  shield  manufactured  from  the  same  material. 

The  structure  of  the  shield  is  shown  as  follows.  In  the  case  of  a  monolithic  shield,  the  thickness  of  the  shield  is  indicated. 
The  notation  2  x  6.0  denotes  the  shield  which  consists  of  two  plates  having  the  same  thicknesses,  6.0  mm.  The  notation 
6.0  +  2.0  denotes  the  two-layer  shield  having  the  front  plate  with  the  thickness  6.0  mm  and  the  rear  plate  with  the  thickness 
2.0  mm.  In  each  column,  the  data  concerning  the  monolithic  shield  is  given  in  the  first  row,  while  data  on  the  layered  shields 
are  shown  below. 

Tables  1-3  contain  data  on  pairs,  impact  velocity-residual  velocity.  Of  two  shields,  the  shield  with  a  smaller  residual 
velocity  and  the  same  impact  velocity  is  assumed  as  having  superior  protective  characteristics.  Tables  4-8  contain  data 
on  the  BLVs  whereby  the  superior  shield  has  a  larger  value  of  the  BLV. 

Comparison  of  the  theoretical  results  obtained  in  Section  3  with  experimental  and  numerical  results  is  shown  in  Tables  9 
and  10.  Inequality  r\res  <  0  or  >jhl  <  0,  where 


Vres  Vres  Vbl  Vbl 

— t. — ,  nb\  =  — — 

V  res  V bl 


(17) 


indicates  the  agreement  between  theoretical  predictions  and  experimental  or  numerical  results  and  ijres  >  0  or  >]hl  >  0  de¬ 
notes  the  discrepancy  from  the  theoretical  predictions.  Inspection  of  Tables  9  and  10  reveals  that  theoretical  predictions 
are  supported  by  experimental  or  numerical  data  in  87%  of  cases. 

Inspection  of  Table  3  shows  that  for  all  values  of  the  impact  velocity,  the  residual  velocity  for  the  two-layered  shield  with 
the  same  thicknesses  of  the  layers  (2  x  4.0  mm)  is  larger  than  for  any  other  partitioning  of  the  shield,  in  compliance  with  the 
theoretical  predictions  in  Section  4  for  the  case  of  convex  function  G(T).  However,  the  BLV  in  Table  8  for  2  x  4.76  mm  shield 
is  larger  than  for  other  two  variants  of  the  shield  with  unequal  thicknesses  of  the  layers.  The  effect  of  the  thicknesses  of  the 
layers  upon  ballistic  properties  of  the  shield  requires  further  investigations.  Indeed,  the  maximum  adverse  effect  of  layering 
with  equal  thicknesses  of  the  layers  in  our  analysis  is  the  consequence  of  the  convexity  of  the  function  G(T),  while  the  avail¬ 
able  experimental  data  afford  also  the  linear  function  G(T),  whereby  the  ballistic  efficiency  of  the  shield  does  not  depend 
upon  the  thicknesses  of  the  layers. 
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Fig.  1.  (a  and  b)  Squared  BLV  vs.  thickness  of  the  shield.  Dotted  lines  are  either  tangents  to  the  curves  at  the  point  T  =  Tmin  or  continuations  of  the  straight 
lines  to  the  left  of  this  point. 


Table  1 

Numerical  results  of  Iqbal  et  al.  [5]  on  residual  velocities  (conical-nosed  projectile,  Weldox  460E  steel  shield). 


No. 

Shield  (sizes  in  mm) 

Impact  velocity  (m/s) 

600 

406 

356 

318 

300 

281 

1 

12.0 

555 

333 

269 

217 

189 

154 

2 

2  x  6.0 

554 

337 

273 

221 

195 

163 

Table  2 

Numerical  results  of  Iqbal  et  al.  [5]  on  residual  velocities  (ogive-nosed  projectile,  1100-H12  aluminum  shield). 


No. 

Shield  (sizes  in  mm) 

Impact  velocity  (m/s) 

150 

113 

97 

83 

82 

73 

66 

57 

1 

1.0 

141 

100 

81 

64 

62 

49 

32 

15 

2 

2  x  0.5 

145 

101 

84 

67 

65 

54 

40 

20 

Comparison  of  different  structures  of  the  shield  7-8  in  Table  4,  6-7  and  1 1-12  in  Table  5  and  4-5-6  in  Table  4  supports  the 
theoretical  predictions  in  Section  5:  increase  of  the  number  of  layers  having  the  same  thickness  impairs  ballistic  properties  of 
the  shield. 
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Table  3 

Experimental  results  of  Almohandes  et  al.  [13]  on  residual  velocities  (7.62  mm  bullet,  mild  steel  shield). 


No. 

Shield  (sizes  in  mm) 

Impact  velocity  (m/s) 

706 

755 

775 

805 

826 

1 

8.0 

446 

489 

528 

611 

647 

2 

2.0  +  6.0 

477 

514 

549 

609 

658 

3 

6.0  +  2.0 

471 

522 

- 

616 

654 

4 

2  x  4.0 

486 

533 

567 

643 

667 

Table  4 

Experimental  results  of  Radin  and  Goldsmith  [14]  on  BLV  (conical-nosed  projectile,  2024-  Aluminum  shield). 


No. 

Shield  (sizes  in  mm) 

BLV  (m/s) 

No. 

Shield  (sizes  in  mm) 

BLV  (m/s) 

No. 

Shield  (sizes  in  mm) 

BLV  (m/s) 

1 

3.2 

95 

3 

4.8 

144 

6 

6.4 

184 

2 

2  x  1.6 

93 

4 

1.6 +  3.2 

130 

7 

2  x  3.2 

160 

5 

3  x  1.6 

124 

8 

4  x  1.6 

158 

Table  5 

Experimental  results  of  Iqbal  et  al.  [6]  on  BLV  (ogive-nosed  projectile,  1100-H12  aluminum  shield). 


No. 

Shield  (sizes 
in  mm) 

BLV 

(m/s) 

No. 

Shield  (sizes 
in  mm) 

BLV 

(m/s) 

No. 

Shield  (sizes 
in  mm) 

BLV 

(m/s) 

No. 

Shield  (sizes 
in  mm) 

BLV 

(m/s) 

No. 

Shield  (sizes 
in  mm) 

BLV 

(m/s) 

1 

1.0 

45.3 

3 

1.5 

54.3 

5 

2.0 

67.1 

8 

2.5 

79.3 

10 

3.0 

90.4 

2 

2  x  0.5 

39.2 

4 

3  x  0.5 

51.0 

6 

2  x  1.0 

62.9 

9 

5  x  0.5 

66.0 

11 

2  x  1.5 

78.6 

7 

4  x  0.5 

61.7 

12 

3  x  1.0 

74.2 

Table  6 

Numerical  results  of  Teng  et  al.  [7]  on  the  BLV  (conical-nosed  projectile,  Weldox  460E  steel  shield). 


Conical-nosed  projectile 

Light  conical-nosed  projectile 

No.  Shield  (sizes  in  mm) 

BLV  (m/s) 

No. 

Shield  (sizes  in  mm) 

BLV  (m/s) 

1  12.0 

305.9 

3 

12.0 

525.9 

2  2  x  6.0 

282.0 

4 

2  x  6.0 

524.5 

Table  7 

Experimental  results  of  Dey  et  al.  [8]  on  the  BLV  (ogive-nosed  projectile,  Weldox  460E  steel  shield). 


No. 

Shield  (sizes  in  mm) 

BLV  (m/s) 

1 

12.0 

318.1 

2 

2  x  6.0 

288.3 

Table  S 

Experimental  results  of  Woodward  et  al.  [15]  (conical-nosed  projectile,  2024-T351  aluminum  alloy  shield). 


No. 

Shield  (sizes  in  mm) 

BLV  (m/s) 

1 

9.53 

462 

2 

3.18  +  6.35 

421 

3 

6.35  +  3.18 

474 

4 

2  x  4.76 

484 

5 

3  x  3.18 

473 

6 

6  x  1.59 

454 

Therefore,  the  available  experimental  and  numerical  results  support  our  theoretical  predictions. 


7.  Concluding  remarks 

On  the  basis  of  general  assumptions  that  are  pertinent  to  penetration  modeling,  we  investigated  the  effect  of  layering  of 
metallic  shields  against  sharp-nosed  rigid  impactors  and  established/explained  a  number  of  general  assertions  concerning 
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Table  9 

Validation  of  theoretical  results  on  the  residual  velocity. 


Table 

Compared  shields 

Vimp  (m/s) 

Vres-Vres 

V,es 

1 

1-2 

600 

+0.002 

1 

1-2 

406 

-0.012 

1 

1-2 

356 

-0.015 

1 

1-2 

318 

-0.018 

1 

1-2 

300 

-0.032 

1 

1-2 

281 

-0.059 

2 

1-2 

150 

-0.028 

2 

1-2 

113 

-0.010 

2 

1-2 

97 

-0.037 

2 

1-2 

83 

-0.047 

2 

1-2 

82 

-0.048 

2 

1-2 

73 

-0.068 

2 

1-2 

66 

-0.250 

2 

1-2 

57 

-0.333 

3 

1-2 

706 

-0.070 

3 

1-3 

706 

-0.056 

3 

1-4 

706 

-0.090 

3 

1-2 

755 

-0.051 

3 

1-3 

755 

-0.067 

3 

1-4 

755 

-0.090 

3 

1-2 

775 

-0.040 

3 

1-4 

775 

-0.074 

3 

1-2 

805 

+0.003 

3 

1-3 

805 

-0.008 

3 

1-4 

805 

-0.052 

3 

1-2 

826 

-0.017 

3 

1-3 

826 

-0.011 

3 

1-4 

826 

-0.031 

Table  10 

Validation  of  theoretical  results  on  the  BLV. 


Table 

Compared  shields 

Vbt-Vbi 

Vb\ 

4 

1-2 

-0.021 

4 

3-4 

-0.097 

4 

3-5 

-0.139 

4 

6-7 

-0.130 

4 

6-8 

-0.165 

5 

1-2 

-0.135 

5 

3-4 

-0.061 

5 

5-6 

-0.063 

5 

5-7 

-0.080 

5 

8-9 

-0.168 

5 

10-11 

-0.131 

5 

10-12 

-0.204 

6 

1-2 

+0.078 

6 

3-4 

+0.003 

7 

1-2 

+0.094 

8 

1-2 

+0.089 

8 

1-3 

-0.026 

8 

1-4 

-0.048 

8 

1-5 

-0.024 

8 

1-6 

+0.017 

ballistic  properties  of  the  shields.  Omitting  the  technical  details,  the  main  conclusions  of  this  study  can  be  formulated  as 
follows: 

(1)  Layering  impairs  the  protective  properties  of  the  shield  against  perforation. 

(2)  Increasing  the  number  of  layers  having  the  same  thicknesses  impairs  protective  properties  of  the  shield. 

Although  our  theoretical  predictions  are  confirmed  by  comparing  with  experimental  or  numerical  results  it  will  be  useful 
to  conduct  a  dedicated  experimental  study  in  order  to  determine  the  range  of  validity  of  our  theoretical  prediction  depend¬ 
ing  on  the  material  properties  of  the  shield,  shape  of  the  impactor,  thicknesses  of  the  layers  and  impact  velocity. 
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Appendix  A.  Proof  of  Claim 


In  order  to  prove  the  claim  we  use  the  method  of  mathematical  induction. 

As  the  first  step,  we  prove  that  the  claim  is  valid  for  the  case  of  two  plates  ( N  =  2). 

Consider  the  function  F(T)  =  TC'(T)  -  G(T).  Since  G(T)  is  a  convex  function,  F'(T)  =  TG"(T)  >  0  and,  consequently, 

TG'(T)  -  G(T)  >  0.  T  >  Tmin.  (A.l) 

Assume  that  7\  ^  T2  (ifTi  <  T2  we  need  only  to  replace  T 1  <->  T2).  Then  using  the  mean  value  theorem  and  Eq.  (A.l)  we 
obtain: 

5  =  G(Tj  +  T2)  -  G(T,)  -  C(T2)  =  [G(Tj  +  T2)  -  G(T,)]  -  G(T2)  =  T2G'(()  -  G(T2)  >  T2C'(£)  -  T2G'(T2 ) 

=  T2[G'(()  -  C'(T2)]  =  T2G"(d),  (A.2) 

where 

h<£<Tl+  T2.  T2<0<£.  (A.3) 

Since  0  >  Tm in,  we  proved  that  <5  >  0. 

At  the  second  step,  we  assume  that  the  claim  is  valid  for  N  plates  (3  =  SN  >  0)  and  prove  that  it  is  valid  for  N  +  1  plates 

(<Wi  >  0).  Using  the  validity  of  the  statement  for  two  plates  with  the  thicknesses  T(1)  +  T<2)  + - b  T(N)  and  T(N+1)  we  can 

write: 


c(£r®j  =  g^t®  +  t(JV+,) J  >c|jf]r®J  +  g(t(N+1)). 

Then  using  the  validity  of  the  claim  for  N  plates  we  obtain: 

(N+l  \  N+ 1  /  N  \  N+l  /  N 

T>j)  -  Y,  G(r0>)  >  G  Z!  r0)  )  +  G(T(N+1))  ~Y,G{T<j))  =G  £  T® 

j=i  /  j=i  \j=i  J  j=i  \j=i 


-^G(T®)  >0. 
j=i 


(A.4) 


(A.5) 


Appendix  B.  The  worst  layering 

In  this  appendix  we  prove  that  the  worst  layering  occurs  when  the  monolithic  shield  is  replaced  by  the  identical  plates 
having  the  same  total  thickness. 

After  eliminating  the  variable  Tm,  the  problem  can  be  formulated  as  follows: 

A(Tm, . . .  ,T(N -V)  =  YjG(T<i))  +  G(Tsum  -  2T0>)  -*■  min  (B-1) 

J=i  j=i 

subjected  to  the  following  constraints: 

£  T®  <  Tsum  -  Tmin,  T®  >  Tmin,  i  =  1 . N  -  1 .  (B.2) 

j=i 


Let  us  prove  that  U  is  a  convex  function.  Formulas  for  first  and  the  second  order  derivatives  of  G  (v,  fi  =  1 , . . . ,  N  -  1 )  read: 


HMSJT=G'(rM) 

-G'(T(N>),  t(N)  =  rsum-^r®, 
j=i 

(B.3) 

u(v,W  _  <9/1  _  ( 

qjWqjW  1 

G"(T(N)) 

G"(T(V))  +  G"(T(n)) 

if 

if  v=  ju 

(B.4) 

Then  the  characteristic  equation, 

H(U)  _  t]  Hd,2) 

fj(  l.N-1) 

HaV  H(2,2)  _ 

W(2,N-1) 

=  0, 

(B.5) 

^(N- 1,2) 

...  H  N  1  v  1  ), 

can  be  written  as  follows: 

J][C"(T<V))-'/]  =  0 
j=i 

(B.6) 
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because  the  last  term  in  Eq.  (B.4)  is  the  same  for  all  and  can  be  dropped  out  in  the  determinant  in  Eq.  (B.5),  and  the 
nonzero  elements  remain  only  at  the  diagonal.  Since  G"(T)  >  0,  all  the  roots  of  Eq.  (B.6)  are  positive  and,  consequently,  A 
is  a  convex  function.  Taking  into  account  that  the  domain  determined  by  the  linear  inequalities  in  Eq.  (B.2)  is  convex,  we 
conclude  that  the  function  A  has  the  unique  minimum.  The  point  of  minimum  can  be  determined  using  the  conditions 
H(v)  =  0(v  =  1 . IV  —  1).  By  virtue  of  Eq.  (B.3)  these  conditions  imply  that 

G'(T|1))  =  G'(T(2))  =  •  •  •  =  G'(Tin^))  =  G'(Tm).  (B.7) 

The  solution  of  Eqs.  (B.7)  is  Tm  =  T(2)  =  •  ■  •  =  T(AM)  =  Tm  =  Tsum/N. 


References 

[1]  Ben-Dor  G,  Dubinsky  A,  Elperin  T.  Ballistic  impact:  recent  advances  in  analytical  modeling  of  plate  penetration  dynamics.  A  review.  Appl  Mech  Rev 
2005;58(6):355-71. 

[2]  Corbett  GG,  Reid  SR,  Johnson  W.  Impact  loading  of  plates  and  shells  by  free-flying  projectiles:  a  review.  Int  J  Impact  Engng  1996;18(2):141-230. 

[3]  Ben-Dor  G,  Dubinsky  A,  Elperin  T.  Applied  high-speed  plate  penetration  dynamics.  Dordrecht  (The  Netherlands):  Springer;  2006. 

[4]  Borvik  T,  Dey  S,  Clausen  AH.  Perforation  resistance  of  five  different  high-strength  steel  plates  subjected  to  small-arms  projectiles.  Int  J  Impact  Engng 
2009;36:948-64. 

[5]  Iqbal  MA,  Chakrabarti  A,  Beniwal  S,  Gupta  NK.  3D  numerical  simulations  of  sharp  nosed  projectile  impact  on  ductile  targets.  Int  J  Impact  Engng 
2010;37(2):185-95. 

[6]  Iqbal  MA,  Gupta  NK.  Ballistic  limit  of  single  and  layered  aluminium  plates.  Strain  2010.  doi:  1 0.1 111/j. 1475-1 305. 2008. 00601.x. 

[7]  Teng  X,  Dey  S,  Borvik  T,  Wierzbicki  T.  Protection  performance  of  double-layered  metal  shield  against  projectile  impact.  J  Mech  Mater  Struct 
2007;2(7):  1309-30. 

[8]  Dey  S,  Borvik  T,  Teng  X,  Wierzbicki  T,  Hopperstad  OS.  On  the  ballistic  resistance  of  double-layered  steel  plates:  an  experimental  and  numerical 
investigation.  Int  J  Solids  Struct  2007;44(20):6701-23. 

[9]  Glazyrin  VP,  Orlov  MYu.  Modeling  of  layered  shields  perforation.  Vychislitel’nyje  Technologii  2003;8(4):143-51  [in  Russian]. 

[10]  Iqbal  MA,  Gupta  NK.  Energy  absorption  characteristics  of  aluminum  plates  subjected  to  projectile  impact.  Lat  Am  Int  J  Solids  Struct  2008;5:259-87. 

[11]  Liaghat  GH,  Moslemi  Naeini  H,  Felli  S.  The  mechanics  of  normal  and  oblique  penetration  of  conical  projectiles  into  multi  layer  metallic  targets.  Iranian  J 
Sci  Technol,  Trans  B,  Engng  2005;29(b2):241-51. 

[12]  Zhou  DW,  Stronge  WJ.  Ballistic  limit  for  oblique  impact  of  thin  sandwich  panels  and  spaced  plates.  Int  J  Impact  Engng  2008;35:1339-54. 

[13]  Almohandes  AA,  Abdel-Kader  MS,  Eleiche  AM.  Experimental  investigation  of  the  ballistic  resistance  of  steel-fiberglass  reinforced  polyester  laminated 
plates.  Compos  Part  B  1996;27(5):447-58. 

[14]  Radin  J,  Goldsmith  W.  Normal  projectile  penetration  and  perforation  of  layered  targets.  Int  J  Impact  Engng  1988;7(2):229-59. 

[15]  Woodward  RL,  Cimpoeru  SJ.  A  study  of  the  perforation  of  aluminium  laminate  targets.  Int  J  Impact  Engng  1998;21(3):1 17—31. 

[16]  Gupta  NK,  Madhu  V.  An  experimental  study  of  normal  and  oblique  impact  of  hard-core  projectile  on  single  and  layered  plates.  Int  J  Impact  Engng 
1997;19(5-6):395-414. 

[17]  Honda  K,  Takamae  G,  Watanabe  T.  On  the  measurement  of  the  resistance  of  shield  plates  to  penetration  by  a  rifle  bullet.  Tohoku  Imperial  University, 
1st  Series  1930;19:703-25. 

[18]  Nixdorff  K.  Discussion  of  two  theories  on  the  penetration  of  multilayer  metallic  targets.  Trans  CSME  1987;  11  (3):  161 -78. 

[19]  Ben-Dor  G,  Dubinsky  A,  Elperin  T.  Ballistic  properties  of  multilayered  concrete  shields.  Nucl  Engng  Des  2009;239(10):1789-94. 

[20]  Recht  RF,  Ipson  TW.  Ballistic  perforation  dynamics.  ASME  J  Appl  Mech  1963;30(3):384-90. 

[21]  Ohte  S,  Yoshizawa  H,  Chiba  N,  Shida  S.  Impact  strength  of  steel  plates  struck  by  projectiles.  Bull  JSME  1982;25(206):1226-31. 

[22]  Gupta  NK,  Ansari  R,  Gupta  SK.  Normal  impact  of  ogive  nosed  projectiles  on  thin  plates.  Int  J  Impact  Engng  2001;25(7):641-60. 

[23]  Showalter  DD,  Placzankis  BE,  Burkins  MS.  Ballistic  performance  testing  of  aluminum  alloy  5059-H131  and  5059-H136  for  armor  applications.  US  Army 
Research  Laboratory.  Report  ARL-TR-4427;  2008. 


